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The coupling of Fe magnetic moments in LaFeAsOi-xFa; with the As A\ g phonon is calculated. 
We present first principles calculations of the atomic and electronic structure of LaFeAsO as a func- 
tion of electron doping. We perform calculations using the virtual crystal approximation as well as 
supercell calculations with F substitutional impurity atoms. The results validate the virtual crystal 
approximation for the electronic structure near the Fermi level. Its is found that the electronic 
density of states at the Fermi level is maximum for x=0.125, enhancing the electron-phonon interac- 
tion. An additional increase of the electron-phonon parameter A is obtained if the coupling between 
the Ai g phonon and the Fe magnetic moment is included. It is found that the electron-phonon 
interaction can be one order of magnitude larger than its value if no spin resolution is included in 
the calculation. The implications of these results on the superconducting transition are discussed 

PACS numbers: 71.20.-b, 71.38.-k, 74.20.-z, 74.70.-b 



Since the discovery of the Fe-based superconductors a 
few years ago [J, Q , several properties common to most 
of the compounds have been well established and general 
agreement has been reached. However, there are very 
important aspects that have not been settled The 
most important one, is the mechanism responsible for 
their superconducting behavior. Early theoretical work 
0, HI ruled out any phonon's role. More recent calcula- 
tions G] reveal that, in the non-magnetic phase, phonons 
coupled with magnetism contribute to superconductivity 
although not enough to explain the high critical tem- 
perature. Another aspect to be clarified involves the Fe 
magnetic moment. The low temperature stripped an- 
tiferromagnetic order in the parent compound is gener- 
ally accepted, although the magnitude of the calculated 
magnetic moment is much larger than the observed one 
• Fluctuations mechanisms as well as orbital magnetic 
ordering 0, 0] have been proposed to explain this dis- 
crepancy. The coupling of the Fe magnetic moments 
with phonons has been clearly established at least in 
CaFe2As2, in the undoped orthorombic phase [10|, and in 
SmFeAsOi-ajFa; ll| . Also the possibility of a non-zero 



Fe magnetic moment in the superconducting phase has 
not been ruled out [12j . 

Here we propose that the magnetic moment in Fe 
atoms does not disappear in the superconducting phase. 
Instead, what disappears is the long-range collinear an- 
tiferromagnetic order while keeping the As-mediated an- 
tiferromagnetic local coupling between the Fe atoms. To 
this end, we have calculated in the manner described 
in [l3| how the energy barrier between two equivalent 
collinear antiferromagnetic configurations varies with the 
extra electron concentration in LaFcAsOi-^Fa,. We as- 
sume that a non-collinear configuration takes place be- 
tween the two minima as discussed in II 31 ] - The results 



indicate that the barrier between two equivalent antifer- 
romagnetic arrangement is very small beyond x=0.1 up 
to the disappearance of the magnetic moment around 
x=0.3. This low barrier can be responsible for the disap- 
pearance of the antiferromagnetic long-range order since 
the system may fluctuate between the two equivalent 
magnetic arrangements. 
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FIG. 1: (Color online) Energy barrier for the transition be- 
tween two equivalent antiferromagnetic stripped collinear ar- 
rangements (left inset) through an intermediate non collinear 
state (right inset) in LaFeAsOi-zFz as a function of x. The 
energy difference per formula unit is shown. 

To study how the electronic structure, and eventually, 
the electron-phonon interaction varies with doping we re- 
strict ourself to the antiferromagnetic collinear arrange- 
ment since the non-collinear one has a similar electronic 



of the calculation are shown in Figure [T] The results 



structure (see [l3j) and its calculation is much more com- 
putational demanding. To obtain the electronic structure 
we have performed density functional }l4lll5l ] calculations 
using the SIESTA code flil [l7| which uses localized or- 
bitals as basis functions 1181 . In our calculation we use 
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a double £ polarized basis set, non-local norm conserv- 
ing pseudopotentials and a local density approximation 
(LDA) for exchange and correlation. The calculations are 
performed with stringent criteria in the electronic struc- 
ture convergence (down to 10~ 5 in the density matrix), 
Brillouin zone sampling (up to 18000 fc-points), real space 
grid (energy cut-off of 500 Ryd) and equilibrium geom- 
etry (residual forces lower than 10 -2 eV/A). Due to the 
rapid variation of the density of states at the Fermi level, 
we used a polynomial smearing method [l9| . To simulate 
the effect of doping we use the virtual crystal approxima- 
tion (VCA) 

The calculated total densities of states for various fluo- 
rine contents are drawn in Figure [5] We hrst observe the 
antiferromagnetic pseudo gap at the Fermi level for x=0. 
As the excess of electrons increases with the fluorine con- 
centration the size of the gap decreases and eventually 
disappears between x=0.25 and x=0.30. It is interesting 
to notice that the antiferromagnetic gap lies fully bellow 
the Fermi level for x > 0.15. Moreover, the peak in the 
density of states at the antiferromagnetic gap upper edge 
crosses the Fermi level between x=0.10 and x=0.15. 



Before we go into the discussion of the electron-phonon 
interaction parameter we analyze to what extent the vir- 
tual crystal approximation is a good starting point to 
study the electronic structure of these materials. To this 
end, we have calculated the equilibrium atomic structure 
and the electronic density of states for various super- 
cells size substituting an oxygen atom by a fluorine one. 
In this manner we have performed "exact" calculations 
for x=0.25, 0.125 and 0.0625. The calculated density 
of states after full relaxation of the atomic positions are 
shown in Figure |3] which shows an excellent agreement 
with the VCA results. The VCA density of states repro- 
duces the main features of the super-cell calculations, in 
particular, the peak in the density of sates crossing the 
Fermi level at around x=0.125. From these results it is 
clear that the details of how impurities are incorporated 
in the oxygen layer does not affect much the electronic 
structure around the Fermi level which is mainly dom- 
inated by the Fe electrons. The main effect of doping 
impurities at the oxygen layer is to pump electrons into 
the As-Fe-As layer, as anticipated [21 1 . 
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FIG. 2: (Color online) Variation of the total density of states 
of LaFeAsOi-^F^ with the F content x. The calculation is 
for the virtual crystal approximation. A small imaginary con- 
tribution has been added to the energy to smooth the curves. 
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FIG. 3: Calculated LaFeAsOi-^Fa; electronic densities of 
states using super-cells ("exact") (panels (a), (b) and (c)) 
and using the virtual crystal approximation (panels (d),(e) 
and (f)) for various values of x. A small imaginary contribu- 
tion has been added to the energy to smooth the curves. 

We have also calculated the influence of an F impurity 
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on the electronic charges and magnetic moments of the 
Fe atoms. In Table Q] we show the results for a supercell 
Fei 6 Asi 6 Lai 6 0i5F 1 which corresponds to x=0.0625. We 
notice that the extra charge is distributed rather evenly 
through the Fe atoms and does not distinguish between 
those closer to the impurity and the other Fe atoms. On 
the other hand, although the magnetic moment of the Fc 
atoms closer to the F impurity display a magnetic mo- 
ment larger than the others, the difference is rather small 
and we can rule out any relevant polaron formation. The 
antiferromagnetic coupling between successive Fe layer 
agrees with the experimental magnetic order of this ma- 
terial. The above results stress the appropriateness of 
the VCA approximation. Also, there is no significant in- 
fluence of the presence of F impurities in the geometry 
of the Fe-As layers next to it. We find that the Fe-As 
distance depends on the extra charge at the As-Fe-As 
layers, rather than on the details concerning the oxygen 
layer and the substitutional impurities. 



TABLE I: Electronic charge and magnetic moment (in Bohr 
magnetons) at the Fe atoms of FeieAsieLaigOisFi in the 
super-cell calculation. The corresponding VCA results are 
indicated. 



Fe Layer 


Qi Q 2 Qs 


Q4 


1 


8.298 8.305 8.312 


8.314 


2 


8.298 8.306 8.300 


8.304* 


3 


8.298 8.305 8.300 


8.304* 


4 


8.298 8.306 8.312 


8.313 


Fe Layer 


Ml M2 M3 




1 


-1.222 -1.219 1.174 


1.212 


2 


1.220 1.218 -1.238 


-1.234* 


3 


-1.222 -1.219 1.240 


1.236* 


4 


1.222 1.220 -1.176 


-1.213 


Qvca- 8.301; fivcA- 1.247 


* Fe atoms closest to the F impurity 



To address the electron-phonon interaction in this 
system and the superconducting A parameter, we have 
considered the symmetric out-of-plane As Ai g mode at 
fen =0, in which the Fe atoms remain fixed whereas 
the As atoms move perpendicularly to the FeAs layers, 
expanding and compressing the Fe-As bonds. For this 
phonon mode the diagonal electron-phonon matrix ele- 
ment can be written [22|: 

< ~t,n\H e - ph {v)\~t,n>= {h/2M a uj u ) 1 / 2 e a „-D(t,n) 
where e° is the polarization vector, u v is the frequency 
of the vibrational mode involved and ~D( k , n) is the de- 
formation potential for the state ( k , n) 

The A parameter associated to a vibrational mode can 
be written as a sum of the contributions of the electronic 
bands crossing the Fermi level in the form: 
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FIG. 4: (Color online) Band structure with frozen Ai g phonon 
(vibration of As atoms perpendicularly to the Fe plane) for 
a doping of 0.10 electrons per formula unit. The results are 
for the virtual crystal approximation. The dashed bands cor- 
respond to zero phonon amplitude and the red (blue) bands 
correspond to an expansion (compression) of the Fa- As planes 
distance u of 0.02 A. (a) Spin independent calculation, (b) 
Stripped antiferromagnetic magnetic order. 



where N n (Ep) is the nth-band density of states at the 
Fermi level, uja 19 the phonon frequency, Ma s the mass of 
the arsenic atoms and D n is the variation of the energy 
band E n at the Fermi energy the phonon with the phonon 
amplitude: 

n _ l dE n (E F ) 
1Jn ~ v/2 du 

Before we calculate the A parameter it is worth look- 
ing at the deformation potential associated to the Ai g 
phonon. Although a supercell approximation could be 
more precise [23[, the calculations are performed within 
the VCA. In Figure |4] we show, for the fluorine concen- 
tration of 10%, how the electronic band structure is per- 
turbed by the presence of the phonon. In the case of a 
non magnetic calculation (Figure H (a)) the effect of the 
phonon in the bands at the Fermi level is small in all the 
bands but one {d xy ). On the contrary, in the case of the 
calculation for the antiferromagnetic configuration (Fig- 
ure |4] (b)), all the bands are perturbed by the presence 
of the phonon, and, in addition, there is a flattening of 
the bands due to the antiferromagnetic pseudo-gap. The 
effect of the A\ g phonon is to modulate the size of the 
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Fe magnetic moment with the vibration amplitude and 
therefore inducing an additional shift to the bands. A 
similar behavior was previously discussed fl3j |. 




magnetic calculation but not large enough to explain by 
itself the superconducting critical temperature. All the 
vibrational modes should be included in the calculation 
to obtain the total A parameter. The connection between 
the Fe-As bond length and the magnetic moments in Fe 
atoms is important (as discussed by Yildirim jiij]) and is 
relevant in the electron-phonon interaction. 

In summary, we propose that in the superconduct- 
ing phase of LaFeAsOi-^F^, the Fe atoms have a finite 
magnetic moment that fluctuates between two equiva- 
lent collinear antiferromagnetic configurations. The pres- 
ence of this moment varies substantially the density of 
states near the Fermi level and enhances dramatically 
the clcctron-phonon interaction at least fort the As A lg 
phonon mode. The maximum of the calculated A param- 
eter takes place at x=0.125. These results suggest that 
electron-phonon interaction, coupled with the Fe mag- 
netic moments, must be carefully revised, before ruling 
out its connection with Fe-based compounds supercon- 
ductivity. 
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I also thank E. Bascones, M. J. Calderon, M.L.Cohen, 
G. Gomez-Santos, S.G. Louie, D. Sanchez-Portal and B. 
Valenzuela for very helpful comments about this work. 
This work was supported by the Spanish Ministry of Sci- 
ence and Innovation through grants FIS2009-12712 and 
CSD2007-00050. 



FIG. 5: (Color on line) Variations of the density of states 
at the Fermi level and parameter A of LaFeAsOi-^Fa, with 
the fluorine concentration x. (a) Density of states for the 
antiferromagnetic collinear order (black squares) and non- 
magnetic (red circles) calculations, (b) Calculated parame- 
ter A for the antiferromagnetic collinear order (black squares) 
and non-magnetic (red circles) . 

We have calculated the A parameter associated to this 
phonon in the approximation discussed above and within 
the VCA. The results, along with the variation of the den- 
sity of states at the Fermi level, are shown in Figure [5j 
where a phonon energy of 25 meV is considered. We ob- 
serve that the antiferromagnetic pseudo-gap induces an 
increase of the density of states at the Fermi energy with 
respect to the non-magnetic calculation (see Figure [5]) . 
This increase would enhance the A parameter by a factor 
of two at the most. However, if the deformation potential 
is included, the A parameter can be one order of magni- 
tude larger than its value for the non-magnetic calcula- 
tion (Figure [5] (b)). In addition, if we consider the renor- 
malization of the parameter A due to the renormalization 

2 

of the phonon mode oj 2 = j-j-§T then A = 1 _^ An which for 
the maximum value at x=0.125 becomes A = 0.16. This 
is indeed a substantial increase with respect to the non- 
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